Optical Kagome lattice for ultra-cold atoms with nearest neighbor interactions 
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We propose a scheme to implement an optical Kagome lattice for ultra-cold atoms with control- 
lable s-wave interactions between nearest neighbor sites and a gauge potential. The atoms occupy 
three different internal atomic levels with electromagnetically-induced coupling between the levels. 
We show that by appropriately shifting the triangular lattice potentials, experienced by atoms in 
different levels, the Kagome lattice can be realized using only two standing waves, generating a 
highly frustrated quantum system for the atoms. 
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Ultra-cold atomic gases in periodic optical lattice po- 
tentials can form a highly-controllable quantum many- 
particle system that has demonstrated interesting analo- 
gies to crystal lattice systems of strongly interacting elec- 
trons. Among such experimental developments are the 
Mott insulator states of both bosonic [1| and fermionic 
atoms (2, Q , fermionic superfluidity Q , atom transport 

Hi Hi ^-iid super-exchange correlations in a collection 
of double wells [1]. The rapid progress has generated 
interest in engineering ultra-cold atomic lattice systems 
that could provide clean realizations of Hubbard mod- 
els of strongly-correlated crystal lattice systems with the 
potential use of atoms as quantum simulators of some 
unresolved models, e.g., in high- Tc superconductivity Q. 

One of the challenging problems in the theory of 
strongly correlated systems has been to characterize the 
phase diagrams of two-dimensional (2D) and 3D Hub- 
bard models of geometrically frustrated lattices, with 
competing interactions resulting in highly degenerate 
ground states. Of particular interest are corner-sharing 
networks of complete graphs, such as Kagome and py- 
rochlore/checkerboard lattices. It was recently proposed 
that frustration in pyrochlore and diamond lattices could 
generate fractional charges in the presence of nearest- 
neighbor (nn) repulsion [10| , with similar excitations po- 
tentially existing in a Kagome lattice. Moreover, spin-1/2 
nn Heisenberg antiferromagnet on a Kagome lattice pro- 
vides arguably the most promising candidate for a quan- 
tum spin liquid, with the Kagome lattice exhibiting a 
hig her degree of frustration than, e.g., a triangular lattice 
[llj l , and it has recently attracted considerable interest in 
Volborthite and Herbertsmithite combounds. Despite ex- 
tensive theoretical effort the true nature of the ground 
state of the Kagome system has been evasive [12] ■ In 
addition, Kagome systems can exhibit, e.g., kinetic fer- 
romagnetism [l3j and production of trimerized and ideal 
Kagome lattices for ultra-cold atoms has also started at- 
tracting theoretical interest [l3| . 

In this paper we propose a constructed 2D optical 
Kagome lattice system of neutral atoms with collisional 
interactions between nn sites. The hopping between the 
nn sites is induced by electromagnetic (em) transitions 
between three internal atomic states. Due to the spin- 
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FIG. 1: Left: In the Kagome lattice the atoms occupy three 
electronic ground states with a common excited state. The 
frequencies of the two lattice lasers are tuned between the 
atomic resonance frequencies in such a way that us > LVa > 
UJ2 > ojp > u)\. The polarization components of the lattice 
lasers that couple to each individual transition are indicated. 
Right: Atoms in each of the three internal levels \j) {j = 
1,2,3) experience an equilateral triangular lattice potential. 
These lattices are shifted with respect to each other, so that 
the lattice site occupied by the atoms in |2) (|3)) is at the 
midpoint of the triangle side of |1) defined by the unit vector 
k^/fc — {By — \/3e2)/2 {bz). Each site has four nn sites, so 
that, e.g., the site of |1) has two nn sites of |2) and two of |3). 



dependent lattice system, the Kagome lattice can be pre- 
pared with only two optical standing waves (SWs), as 
compared to six SWs in previous proposals |14{ . with 
the additional advantage of controllable non-local two- 
body interactions between different sites and the possi- 
bility for the creation of Abelian and non-Abelian guage 
potentials. The strength of the em-driven hopping can 
be tuned over a wide range of values with respect to 
tunneling between more distant sites and the nn colli- 
sions. Spin-dependent lattices where the different atomic 
spin components were moved around inde pen dently was 
experimentally created using ®^Rb atoms jl5|. Alkaline- 
earth-metal atoms and rare-earth metals with narrow op- 
tical resonances (e.g. Sr, Yb) are particularly suitable for 
realizing spin-dependent optical lattices, because of slow 
loss rates due to spontaneous emission [T^ . 

We consider ultra-cold (bosonic or fermionic) atoms 
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occupying three internal sublevels |j) (j = 1,2,3) of the 
same atom that are coupled by em transitions. The atom 
dynamics is assumed to be restricted to a 2D layer in the 
yz plane due to a tightly confined 2D pancake-shaped 
trap. Such a confinement may be experimentally cre- 
ated, e.g., using magnetic or optical fields Along 
the 2D trap each species experiences a triangular optical 
lattice potential, generated by two SWs with wavcvcc- 
tors kjj.^ = k{ey ± •\/3e2)/2. The atom-lattice model is 
based on an atomic (single-band) Hubbard Hamiltonian 
[l9t where the atoms occupy the lowest mode of each 
lattice site. We show that by appropriately tuning the 
frequencies of the lattice lasers, the lattice potentials of 
the three species can be shifted in a triangular shape to 
form a Kagome lattice pattern. Then the four nn sites are 
occupied by the atoms in the other two sublevels and the 
hopping of the atoms between adjacent sites, with ampli- 
tude Kjfc, only occurs as a result of driving by coherent 
em fields that change the internal atomic level. The hop- 
ping to the nearest sites occupied by the atoms in the 
same sublevel results from tunneling between the sites 
and, in sufficiently deep lattices, it may be suppressed. 
Moreover, the specialty of the proposed scheme is that 
by adjusting the overlap of the nn lattice site wave func- 
tions (Wannier functions) we can prepare a lattice system 
with a non-negligible, controllable s-wave interaction be- 
tween adjacent sites [j ^ k), with the both limits 
Ujj^ ^ Kjk and J7™ ^ Kjk achievable, providing a frus- 
trated quantum system with nn interactions. 

As a realization of a frustrated Kagome lattice with de- 
sired interactions, we consider a tripod four-level scheme; 
Fig. [T] The atoms in the three electronic ground states 
are coupled to a common electronically excited level |e) 
by different resonant transition frequencies ujj {ujs > 
LU2 > LOi). The lattice laser a, with the frequency uja, 
is blue-detuned with respect to two of the transitions 
(5q,j = ujj — uJa < 0, for j = 1,2, and red-detuned with 
respect to one 6 a, 3 > 0. Then the atoms in the states 
|1), 1 2) (1 3)) arc attracted towards low-intensity (high- 
intensity) regions of a. The lattice laser /3, with the 
frequency cjp^ is blue-detuned from one of the transitions 
(5/3^1 < and red-detuned from the others. For simplic- 
ity, we assume |e) = |e, m), |1) = |1, m — 1), |2) = |2, m), 
and |3) = |3, m + 1), so that the transitions |1) |e), 
1 3) |e), and |2) |e) have dipole matrix elements 
del = (e|d|l) = S)(e|l,l)e;i, d^a = S)(e| - l,3)eli, 
and de2 = S(e|0, 2)eQ, coupling to the light with po- 
larizations (7+, <7~ , and eg, respectively. Here J) is the 
reduced dipole matrix element. {e\ag) arc the Clebsch- 
Gordan coefficients, and we use e+i = — (e^. -|- iey)/^/2, 
e_i = (e^i, — iey)/^/2, and eo = e^. We assume the po- 
larizations of the two SWs to be orthogonal • = 0, 

so that the lattice potentials read Vj ~ V°' + with 

V," ^Sj^Er sin^[k{y + V3z)/2 -t- ipj], 
Vf ^s.pEr sm^[k{y - ^/iz)/2 + 77,] , (1) 
Vi = = <<23 = 0, <^2 = m = % = 7r/2. Here s^^ oc 



|(de_,- • e^)^/5^j| [t] = a,/3) denotes the lattice strength 
in the lattice photon recoil energy units Ej. = /i^fc^/2m, 
depending on the light polarization, atomic sublevel, and 
the detuning. To produce three triangular lattices the 
lasers need to couple simultaneously to all the three tran- 
sitions, so that Sja and sjp are non- vanishing for all j. 
This can be obtained by choosing the polarization vec- 
tors for the SWs so that • e^.]^ q 7^ together with 
e,, • = (r/ = a, /3) and Gq, • eg = 0. 

We will next demonstrate an example to show that 
such a solution can be found for the polarization vec- 
tors and how selecting the lattice light polarization can 
also be used to control the relative strengths of the three 
triangular lattices. This can be especially useful if the 
absolute values of the detunings are very different. We 
consider the level scheme shown in Fig. [T] and, for sim- 
plicity, assume that \da,3\ = l^aM = 3\6a,i\ = ^ISp^sl = 
\Sf3.2\ = \S[3,i\ (indicating that lui — UJ2 = UJ2 — (^3 and 
that the lasers are detuned exactly at the midpoint be- 
tween the nearest transitions) and that the Clebsch- 
Gordan coefficients are equal. The polarization vectors 
(orthogonal to the corresponding wavevectors ) for the 
SWs are = Qai^VSey + ez)/2 + bae^ and e^j = 
O-piV^By + ez)/2 + bfjBx- Here we choose the complex 
coefficients aa^f3 and ba.fs in such a way that ec, • = 
and \e*_i ■e^l'^ = 3|e^;^ -ed^ and je^^ -e/jp = 3|el;^ -e/jp. 
The last two conditions ensure that we have sir/ = 53^, 
compensating for the different values of the detunings. 
A straightforward algebra yields solutions for aa,f3 and 
ba,(3 with S2ri = 4s3,,/5 (77 = a,/3), resulting in only very 
small differences between the lattice strengths. We also 
find that the two SWs in Eq. ([T|) have equal amplitudes 
for each level, i.e., Sja = Sjfs for all j. 

For each species the lattice potential Vj, in Eq. ([T]) is 
an equilateral triangle with the side of length 2d, where 
d = Tr/^/3k denotes the nn separation; Fig. [T] The min- 
ima of the potential are at the corners of the triangles 
which for |1) are at (y, z) ~ [{n + m)7r/fc, (n — m)d], for 
|2) at [{n + m + l/2)TT/k, {n - m - l/2)d\, and for |3) at 
[{n + 171 + l)7r/fc, (n — m)d], with n, m integers. The tri- 
angular lattices of |2) and |3) are shifted to coincide with 
the side midpoints of the triangular lattice of |1). The 
combined system of interlaced triangular lattices forms a 
Kagome lattice, so that each site has four nn sites which 
are occupied by the two other species (two of each) . 

In order to estimate the relative strengths of the dif- 
ferent terms in the lattice Hamiltonian we evaluate the 
corresponding integral representations. The direct tun- 
neling amplitude, where atoms remain in the same hy- 
perfine level during the hopping process, reads ~ 

-Jdydz (-|^0LVy<Afe«' +'^L^b'/'h"') > 0, where V| = 
dy + d1 and Vf, is given by Eq. H]). The Wannier func- 
tions for atoms in \h) at site n are 4>bniy^ z) and may be 
approximated by the ground state harmonic oscillator 
wave function with the trap frequencies Uy ~ \/2sEr/h 
and Uz — \/6sEr/h which are obtained by expanding 
the optical potential at the lattice site minimum [l9j . 
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The site n' here refers to the nearest site to n occu- 
pied by the same atomic hyperfine level along the di- 
rection of p. where p takes the values of the three tri- 
angle sides: k^^^jZ. Here = J^'' by symmetry. For 

° and J§ we have 4>b,n'iy, z) = 4>b,n{y + VSd, z + d) 
and 4>ii,n'{y,z) = (pb,n{y,z + 2d), respectively. Due to 
the anisotropy of the individual lattice site wavefunc- 
tions, the hopping amplitudes along the z direction differ 
slightly from those along the direction of the two SWs. 
Although a more rigorous calculation of the hopping am- 
plitudes would involve a full band-structure calculation, 
here it is sufficient to provide order-of-magnitude analytic 
estimates using the Gaussian approximations to t/)h„. 

The em field changes the internal level of the atom. 
Because the lattice site minima of the different sublevels 
are shifted with respect to each other, the atoms simul- 
taneously also undergo spatial hopping along the lattice. 
The hopping amplitude reads 

K&c = Jd^r<pl^hTZbc(l)cn', (2) 

where TZbc denotes the effective Rabi frequency for the 
transition between the levels |6) and |c) and n' refers to 
the nn site of species |c) to the site n occupied by the 
species with Kbc proportional to the spatial overlap 
between the atoms in the nn sites. As shown in Fig. [U 
ni2, K23, and K13 represent hopping along the directions 
of k/3, ktj, and z, respectively. As for the direct tunneling 
we then have K12 = K23, but K13 is not exactly equal. 
Note that Kbc in Eq. ^ , unlike the direct tunneling, can 
take positive, negative, or even complex values. 

The explicit expression for the Rabi frequency TZbc de- 
pends on the particular form of the em coupling between 
the internal levels which can be one or multi-photon tran- 
sition. For a two-photon transition via an off-resonant 
intermediate level |e') (which for a laser can be electroni- 
cally excited and for a rf /microwave an electronic ground 
state) we may adiabatically eliminate |e') [2^. We then 
obtain TZbc = Sb^c'^e'bdlic/i'^h'^A) and a contribution 
to the em-induced level shifts ~\£j\'^\dc'j\'^/{2h?A) for 
Here dg'j = ^e'j • and we have assumed that |6) 
is coupled to |e') by the em field with the positive fre- 
quency component Ej^ = ^eb£be^^''''^e'^^''* and detuning 
A. Using the Gaussian approximation to (/)fc„, we obtain 
^ hTZjkejk exp [-(3 + V3)n^y/sJ^/'i8V2], for (j, k) = 
(1,2), (2, 3) and K13 ~ ;l7^l3el3 cxp (-tt^ ^11^/4^6), 
where Sjk = '^SjSk/i^/sJ + a/s^)^ and e^k = 
i^'jk/i^j^i^)]^^'^- Here we have assumed that the lattice 
potentials by the lasers a and /3 in Eq. ([T]) have the same 
amplitude in each level i.e. Sja = Sjp, and we have 
suppressed in the notation the index referring to the par- 
ticular laser. 

The onsite interaction term may also be obtained from 
the Wannier functions Ujj ~ g^^ J dydz \ 4>jn\'^- The 2D 
nonlinearity / mix V27r is given in terms 

of the scattering length Ujk and the 2D trap confine- 
ment Ix = ^/h/mtOx where the oscillator frequency per- 



pendicular to the lattice is ujx- The additional density- 
dependent contribution to the scattering length in 2D 
is negligible when ^/2Trlx/ajj S> \n.{8Tr^^^lxn2DO,jj) [2l[, 
where n2D denotes the 2D atom density. For fermionic 
atoms Qjj = and the onsite interaction term vanishes 
for a single-species gas if the atoms only occupy the 
lowest mode of each lattice site, but is non-vanishing 
for a two-species gas trapped in each lattice site. The 
nn interaction due to the s-wave scattering between the 
atoms in different sublevels is always non-vanishing for 
both bosonic and fermionic atoms. It is generated in 
the spatial overlapping area of the adjacent sites and 
may be estimated by J7"^ ~ 252^^ / dydz \4>jn\'^\4>kn'\^ ■ 
We obtain Ujj ~ ajjS^/'^Er^/lxV^, C/™ ^ 
2ai33'^/'^Ery/s^e-''^^^^^/lxVTr, and for other states 
f/j^" ~ 2a,,3i/4iJ^^e-(3+v^)-'v^/24V^/;,V^. 
The Hamiltonian for the atomic system then reads 

H = Y. ["^^l^" + UkkclclckCk] - (^M^^ + H-c.) 

k (jk)' 

~ ^ [{Kbc^c, + H.c.) + U^^c]c,cick] , (3) 

where (jk) denotes the summation over adjacent lattice 
sites and {jk)' over the nearest sites that are occupied 
by the same atomic species. The level shifts and the 
detunings due to cm-induced hoppings are included in 
Efe. Here > Q but all other coefhcients can take either 
positive or negative values, and Kbc can be complex. 

We may compare the different terms in the Hamilto- 
nian. For typical experimental parameters for bosonic 
atoms we have Ujj ^ C/™ (j k). For simplicity, as- 
suming si = S2, we obtain for lattice heights s = 25 
and 40, U^2/U3j ^ lO'^au/ajj and IQ-^au/ajj. At 
the same lattice heights the em-driven hopping ampli- 
tudes in terms of the direct tunneling and the nn in- 
teractions arc given by ki2/J^°' ~ bAQ{hTZi2 / Er) and 
5100(;^7^l2/S.), and Ki2lU^i ~ 4.2{lx/ai2){hTZi2/ E,.) 
and 8.3{lx / ai2){hTZi2 / Er) ■ In shallow lattices, with weak 
em-coupling between the sublevels hTZi2/Er ^ 1, the nn 
hopping Kbc and the direct tunneling between more 
distant sites may be comparable. In deeper lattices and 
for stronger Kbc the direct tunneling terms may be ig- 
nored. If the transverse confinement of the 2D lattice 
is weak {Ix large) and the Rabi field sufficiently strong 
KRi2lxl{0'i2Er) » 1, we may also neglect C/™ with 
Kbc ^ and in Eq. Q we only keep terms pro- 

portional to Efc, Ukk^ and Kbc- An especially interest- 
ing property of the proposed scheme is that we may 
also find a wide range of parameter values for which 
^bc ^ ^bc- This limit can always be achieved with suf- 
ficiently weak em coupling. If we also simultaneously 
require that Kbc 3> , we may, e.g., at s = 50 se- 
lect hTZi2/Er ~ 5 X 10~'* and Ix — 17ai2, resulting in 
t/™ ~ lOKfcc- For instance, for ®'^Rb the s-wave scatter- 
ing length between \F = 1, Mp = — 1) and |2, -|-1) hyper- 
fine states is about 5.191nm corresponding to the 
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transverse trap frequency of ~ 27r x 15 kHz, achiev- 
able, e.g., by an optical lattice. The interspecies scat- 
tering length could be increased by Feshbach resonances, 
further enhancing the effect of the nn interactions J7™, 
so that one could reach Ujj » J7™ ^ Kbc ^ also in 
shallow lattices. Interspecies Feshbach resonances were 
observed also in ^^Rb between different hyperfine states 
[23| and in optical lattices at low occupation numbers 
the harmful three-body losses are suppressed. More- 
over, alkaline-earth-metal and rare-earth-metal atoms 
with very narrow optical resonances fl^ may allow the 
lattice lasers to be tuned close to the atomic resonance 
making it easier to produce deep lattices. 

The simplest situation is to select the phases of the 
em fields in Eq. ^ so that all the hopping amplitudes 
Kf,c are real and positive. It is, however, also possible 
to engineer a non-uniform phase profile for the hopping 
amplitudes which was in Ref. [2^ 1 proposed as a mecha- 
nism to construct topologically non-trivial ground states 
with fractional fermion numbers in ID. In a 2D lattice 
the technique can be used to create an effective mag- 
netic field for neutral atoms [1^. Here we may simi- 
larly induce a phase for atoms hopping around a closed 
path in the lattice, mimicing a magnetic flux experienced 
by charged particles. We write Kbc = |Kfec|e"^'"^ where 
the phases i'6c(r) may be constant or spatially varying 
[U, |2^ . The hopping around one unit triangle then gen- 
erates the phase Av = 1^12 -I- 1^23 ~ ^^13 for the atoms, 
corresponding to a magnetic flux $ oc Av through the 
area enclosed by the triangle. For instance, for spatially 
constant i>bc with Av ^ 0, the entire lattice area may 
be divided into side-sharing triangles where each adja- 
cent triangle experiences the flux with the opposite sign. 
In the case of atoms occupying more than one sublevel 
in each lattice site, we may also generate non-Abelian 



vector potentials similarly to Ref. [2^ ; see also [13, HI] . 

For strong nn interactions with 1/3 filling, even with- 
out a vector potential, our model produces a frustrated 
ground state where one atom in each triangle of sites 
is strongly influenced by the atoms in other corner- 
sharing triangles. It is helpful to consider a honey- 
comb lattice, formed by connecting the centers of tri- 
angles of the Kagome lattice, where the sites are cou- 
pled by ring- hopping processes [2^. A lattice system 
described by an analogous quantum dimer model on a 
pyrochlore/checkerboard lattice and on a 3D diamond 
lattice was recently shown to support fractional charges 
[lo| . Similar fractional excitations are expected to ex- 
ist in the Kagome lattice system, where they can act 
as independent, deconfined particles over finite distances 
at temperatures above the ordering transition driven by 
quantum fluctuations - in this case by the ring-exchange 
processes [sO] . Atomic states in the prepared lattice sys- 
tem could potentially be detected optically [3l| . 

Our formalism considers one atomic species per lat- 
tice site. It is straightforward to generalize it to the sit- 
uation where a two-species gas is trapped in each site 
with em field inducing hopping for both species. Then 
the onsite interaction and the hopping terms can be ex- 
pressed as an effective Hciscnberg spin-1/2 Hamiltonian 
Hctf ^ Ei^,J)[tzS^S^ ± t^iSfSJ + SfS^)] where the + 

(— ) sign refers to fermionic (bosonic) atoms and 5*^ de- 
note the spin matrices (3^ . The fermionic version has 
been extensively studied in 2D Kagome lattices where 
the ground state of the SU(2) symmetric case {t^ = t±) 
still has unsettled questions, e.g., in the existence of spon- 
taneously broken symmetries and finite energy gaps [12| . 

We are grateful to N. Shannon for explanations of the 
importance of Kagome lattices and EPSRC for funding. 
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